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New shell model Hamiltonians are derived for the T = 1 part of the residual interaction in the 
f5/2 P3/2 Pi/2 g9/2 Hiodel space based on the analysis and fit of the available experimental data for 
28Ni29-2iNi5o isotopes and 29CU50- so^Snso isotones. The fit procedure, properties of the determined 
effective interaction as well as new results for valence-mirror symmetry and seniority isomers for 



nuclei near ^*Ni and ^""Sn are discussed 
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Neutron-rich nickel isotopes in the vicinity of 28Ni5o 
are currently in the focus of modern nuclear physics and 
astrophysics studies 0, SS 0, S S . The enormous in- 
terest in this region is motivated by several factors. The 
primary issue concerns the doubly magic nature of 28-'^i 
and understanding the way in which the neutron excess 
will affect the properties of nearby nuclei and the ^*Ni 
core itself. The shell-model orbitals for neutrons in nu- 
clei with Z = 28 and N=28-50 (^^Ni-'^^Nj) ^-^e same 
as those for protons in nuclei with N=50 and Z=28-50 
('^*Ni-i'^°Sn). Thus it is of interest to understand the 
similarities and differences in the properties of these nu- 
clei with valence- mirror symmetry (VMS) The as- 
trophysical importance is related to the understanding 
of the nuclear mechanism of the rapid capture of neu- 
trons by seed nuclei through the r-process. The path of 
this reaction network is expected in neutron-rich nuclei 
for which there is little experimental data, and the precise 
trajectory is dictated by the details of the shell structure 
far from stability. 

Experimental investigations of neutron-rich nuclei have 
greatly advanced the last decade providing access to 
many new regions of the nuclear chart. Nuclear struc- 
ture theory in the framework of shell-model configuration 
mixing has also advanced from, for example, the elucida- 
tion of the properties of the sd-shell nuclei (A = 16 — 40) 
in the 1980's ^ to those of the pf^shell {A = AQ - 60) 
in current investigations 0, . Full configuration 
mixing in the next oscillator shell (sdg) is presently at 
the edge of computational feasibility. For heavy nuclei 
the spin-orbit interaction pushes the gg /2 and /y /2 orbits 
down relative to the lower-Z orbits. Thus the most impor- 
tant orbitals for neutrons in the region of ^^Ni to ^^Ni are 
P3/21 /5/21 Pi/2 B.nd g^ji (referred to from now on as the 
p/5 /259 /2 model space) . It is noteworthy that this model 
space is not affected by center-of-mass spurious compo- 
nents. Full configuration mixing calculations for neutrons 
or protons in this model space are relatively easy. The 
work we describe here on the T = 1 effective interactions 
will provide a part of the input for the larger model space 
of both protons and neutrons in these orbits where the 



maximum m-scheme dimension is 13,143,642,988. This 
proton-neutron model space is computationally feasible 
with conventional matrix-diagonalization techniques for 
many nuclei in the mass region A=56-1GG, and Quantum 
Monte Carlo Diagonalization techniques ^13] or Exponen- 
tial Convergence Methods 0| can be used for all nuclei. 

The present paper reports on new effective interactions 
for the p/5/239/2 model space derived from a fit to exper- 
imental data for Ni isotopes from A = 57 to A = 78 and 
TV = 50 isotones from ™Cu to ^^^Sn for neutrons and pro- 
tons, respectively. Predictions for the ^^~^^Ni isotopes 
are made using the new effective interaction. For the first 
time the calculated structures of the 68,70,72,74,76]\jj jgQ_ 
topes and the ^"Zr, ^■^Mo, ^''Ru, ^^Pd, ^^Cd are compared 
and analyzed with respect to the VMS concept jS] . Our 
work provides a much improved Hamiltonian for Z = 28 
over those considered in smaller model spaces flU^lsHTsj. 
and also provides a new Hamiltonian for iV = 50 that is 
similar to those obtained previously (3 ^3 ■ 

The effective interaction is specified uniquely in terms 
of interaction parameters consisting of four single-particle 
energies and 65 T = 1 two-body matrix elements 
(TBME). The starting point for the fitting procedure 
was a realistic G-matrix interaction based on the Bonn-C 
NN potential together with core-polarization corrections 
based on a ^^Ni core j^^. The low-energy levels known 
experimentally are not sensitive to all of these parame- 
ters, and thus not all of them can be well determined by 
the selected set of the energy levels. Instead, they are 
sensitive to certain linear combinations of the parame- 
ters. The weights and the number of the most important 
combinations can be found with the Linear Combination 
Method (LCD) [l|. Applying LCD for our fit we found 
that convergence of the in the first iteration is achieved 
already at 20 linear combinations and we have chosen this 
as a reasonable number for all following iterations. We 
performed iterations (about six) until the eigen-energies 
converged. The values of the neutron interaction param- 
eters are adjusted to fit 15 experimental binding energies 
for ^''"''^Ni and 91 energy levels for ^""''■^Ni. The nu- 
clei below ^"Ni were not emphasized in the fit due to 
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FIG. 1: Calculated and experimental excitation energies of 
the 2f and 4^ states in A = 58 — 76 even-even nickel isotopes. 
Calculated levels are given by circles connected by the dashed 
line. Experimental data are depicted by squares connected by 
the solid line. 



the increased role of excitations from the /7/2 orbit as 
^^Ni is approached 12]. In the absence of experimental 
data on the binding energy of nuclei near ^^Ni we include 
the SKX 20] Hartree-Fock value of -542.32 MeV for the 
binding energy of ^^Ni as a "data" for the fit. Our cal- 
culated binding energies for ^^~^^Ni isotopes agree well 
with the recent corresponding extrapolations from Ref. 
|2ljl . For protons 19 binding energies and 113 energy lev- 
els were used in a fit (this data set is similar to that used 
in Ref. The average deviation in binding and ex- 

citation energies between experiment and theory is 241 
keV and 124 keV for neutrons and protons, respectively. 
A detailed report of the new interactions will be given 
elsewhere ^22] . 

In this letter we emphasize some of the interesting re- 
sults for known nuclei and the extrapolation to properties 
of unknown nuclei. To illustrate some general properties 
of the new interactions we plot the excitation energies 
of the 2i and the 4^ states for neutrons and protons 
in Figs. 1 and 2, respectively. The systematics shows 
good agreement between shell-model calculation and ex- 
periment. There is some similarity in the trends for the 
nuclei with A = 68 — 76 and A = 90 — 98, that is referred 
to as the VMS However, the left part of Figs. 1 and 
2 (yl = 58 - 66 for nickel isotopes and A = 80 - 88 for 
= 50 isotones) are drastically different. Two nuclei, 
^^Ni and *^Sr, show the most profound differences in the 
location of the 4^ state. The energy gaps between the 
4^ and the 2f states in the ^^Ni and the ^''Se are also 
obviously distinct. 

These differences may be qualitatively understood 
from the ordering of the single-particle energies (SPE) 
for both cases (see Fig. For neutrons the lowest 
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FIG. 2: The same as in the Fig. [T]for A = 80 - 
A'' = 50 isotones. 
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FIG. 3: The neutron ('^''Ni) and proton (^''Cu) single-particle 
energies (SPE) relative to the ^''Ni and the ^**Ni cores, respec- 
tively. The SPE for neutron holes in ^*Ni and proton holes 
in ®^In are also shown. The SPE values relative to the cor- 
responding cores (p3/2 orbital) are given below (above) the 
plotted lines. The relative SPE for ^^Cu and ^®Sn are similar 
to those of the mirror nuclei ^'^Ni and ®^In, respectively. 



orbital is P3/2, which is followed by the /5/2, P1/2 f^nd 
gg/2 orbitals. This ordering is similar to the familiar 
cases of interactions in the p/-shell. For the protons 
we obtain the f^^2 orbital as the lowest similar to the 
previous Ji and Wildenthal interaction One notes 

that the spacing between P3/2, P1/2 59/2 is rather 
similar in both cases. The fact that the orbital is 
pushed down in energy in ^^Cu relative to ^''Cu may be 
attributed to the neutron mean field of the ''^Ni core. 
The strongly binding monopole interaction in the proton- 
neutron (ttu) spin- flip configuration 7r/5/2i'<?9/2 as com- 
pared to 7rp3/2J^39/2 causes a dramatic down sloping of 
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FIG. 4: Yrast level schemes for even-even Ni isotopes with 
A = 68 — 76. Calculated levels on left side and experimental 
levels on the right side. 



the 7r/5/2 level in Z = 29 (Cu) isotopes upon filling of 
the orbit [2ll23|. 

The difference in the ordering of the proton orbitals as 
compared to neutrons is the main reason for the differ- 
ences (e.g. for the 4+ states) observed in Figs. 1 and 
2. The SPE's impact ground states as well: the /I/2P3/2 
component (/5/2 andp3/2 are filled) constitutes 59.8% for 
the ®*Sr and only 21.4% for the ^^Ki (the VMS partner 
of the *^Sr). This difference determines what happens 
beyond the ^^^Ni ( e.g. ^^m, see also Refs. HQ) or ^SSr 
upon filling the pi/2 and the orbitals. 

To compare the details of the low energy spectra of 
the even 68-76]v^j isotopes and even A 90 - 98 TV = 50 
isotones we show the calculated and experimental ener- 
gies for some levels of interest in Figs. 0]and|Sl respec- 
tively. The energies of the 2^ , 4j^, 6^ and 8^ states are, 
approximately, the same in all four '^""'^^Ni nuclei. A 
similar situation holds for the four isotones ^^Mo-^^Cd, 
where the validity of the generalized seniority approxi- 
mation is well established l25|. Indeed, the calculated 
structure of the wave functions indicate a large contri- 
bution of the gg/2 orbital for these nuclei. However, this 
contribution is not so large to conclude the dominance for 
all four nuclei. This is especially well illustrated by the 
structure of the ground states: the [gQ/2]Q+^^ component 
in the 0^ wave functions for 70,72,74,76]s^j jg 44 %^ 53%^ 
67% and 83%, respectively. It is obvious that, in contrast 
to the single 519/2 orbital approximation, the structures 
of ™Ni and ^^Ni are significantly different. The con- 
tributions of the [39/2] ^-T^* component to the 2f states 
have approximately the same weight. The other compo- 
nents of the wave functions play a very important role. 
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FIG. 5: The same as in the Fig. 21 for even-even A'' = 50 
isotones with A = 90 - 98. 



For instance, the difference between the effective neu- 
tron 52^2; = and gl/^; J = 2 TBME's is 0.373 MeV, 
however the 2+ — 0"*" energy gap in ™Ni or in ''^Ni is 
1.2-1.1 MeV. Thus the largest contribution to the gap 
(0.8-0.7 MeV) is mixing with other configurations. The 
gQ/2 wave function content of the corresponding valence 
mirror partners ^^Mo-^*Cd is slightly larger : 51 %, 60%, 
71% and 84% of [59/2]^-^^" in the ground states of each 
A-isotone, respectively, but the overall situation is rather 

similar to 70,72,74,76nj^ 

The energies of the 2f states for the nickel isotopes are 
systematically lower (~ 0.3 MeV) than for the iV = 50 
isotones. This effect is due to the properties of effective 
(7g^2 TBME's, since the energy gap between the effective 
gl/^; J = and g^^^; J = 2 TBME's is 0.727 MeV and 
0.373 MeV for protons and neutrons, respectively ( i.e. 
they differ by approximately the same ~ 0.3 MeV) and 
the calculated structure of the ^-'^^jv^j g^^^ 92]yjQ_98(^^ 
is similar (the energy gap in the starting renormalized 
Bonn-C Hamiltonian is 0.538 MeV). We interpret this as 
an indication that the Z = 28 proton shell gap near ^*Ni 
is relatively weak compared to the = 50 neutron shell 
gap near ®*Sr. Thus, the nuclei ^^"''^Ni have substantial 
amounts of proton core excitations that are not included 
explicitly in the model space, but are implicitly taken 
into account by the effective TBME. One would like to 
treat the nickel isotopes in a model space which explicitly 
includes the proton excitations, but this is presently not 
feasible in terms of computational power. 

The lowering of the effective J ~ 2- J — gap for the 
nickel isotopes compared to A = 50 has important con- 
sequences for non-yrast states. This is illustrated by the 
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TABLE I: Calculated B(E2; J[ JJ) values for the A = 
70 — 76 Ni isotopes. A reasonable value of 1.0 is assigned, ten- 
tatively, to an effective quadrupole charge en ■ B(E2) values 
are given in units of e^fm*. Calculated and available experi- 
mental lifetimes r for the 8^ state are given in last two rows. 
Experimental excitation energies were used in lifetime calcu- 
lations for ^°Ni . For ^^'''^'^^Ni isotopes theoretical excitation 
energies were used. 
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values for A — 92 — 98, A'' = 50 isotones. A reasonable value of 
2.0 is assigned, tentatively, to an effective quadrupole charge 
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properties of the 6+ states. For ^"^Ru and ^^Pd the sec- 
ond 6"*" state is dominantly seniority v = 4 and Ues above 
the 8^ state, while in ^^Ni and '^''Ni it is well below the 
8i state and is almost degenerate with the 6^ seniority 
1/ = 2 state. Despite the very small splitting between two 
6^ states with dominant seniority 1^ — 2 and v — A they 
are only slightly mixed. The structure and location of 



the states has important implications for the isomeric 
properties of the 8+ states. The B(E2; 8^ — > 6^2) values 
are given m Table |T| The E2 transition between the 8+ 
and 6+ states of the same seniority is forbidden in the 
middle of the shell. It is well known that this senior- 
ity selection rule leads to the isomerism of 8^ states in 
^'*Ru and ^^Pd ^2^, with measured lifetimes r of 102(6) 
fis and 3.2(4) /xs, respectively. Our calculations result in 
lifetimes of the order of fj,s as well. The discrepancy be- 
tween theoretical and experimental B(E2) values for the 
^^Ru is relativel y la rge, and this is common to previous 
calculations [T7Ll25t . 

Turning back to Ni-isotopes, one notes, that pushing 
down of the 6^ = 4 states opens up a new channel for 
the fast E2 decay of the 8+ states that results in the dis- 
appearance of isomeric states in ^^'^*Ni -the correspond- 
ing lifetimes move down to the ns region (see Table HJ. 
Our results with the newly derived interactions fully sup- 
port the above explanation of the absence of an isomeric 
8+ state in ^^Ni proposed in |3|,|^ on the basis of single 
gg/2 orbital shell- model calculations since our calculated 
wave functions of the 8"*", 6+ states in ^^'^^Ni are domi- 
nated by the g^^^ (70%) and g^^^ (85%) configurations, 
respectively. However, it is shown above that this is not 
the case for the and 2^ states. 

Clearly more detailed experimental studies of the nu- 
clei in the vicinity of the ^*Ni are required to verify the 
shell-model predictions. The next step will be to com- 
bine these new neutron-neutron and proton-proton effec- 
tive Hamiltonians with a proton-neutron Hamiltonian to 
describe the wide variety of spherical, collective and co- 
existing features for the ^ = 56 — 100 mass region as 
well as to apply the wavefunctions to the calculations of 
weak-interaction and astrophysical phenomena. 
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